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First, we consider galaxy formation from the viewpoint of hierarchical clustering theory 
i— H ' and discuss the possibility that inhomogeneous models with a local void may be compatible 

with the observed homogeneity of galactic distributions found in recent redshift surveys, 

because their inhomogeneity can be weakened by the difference in the feedback system of 
t-H ' galaxy formation between the inner and outer regions. Next, it is shown with the results of 

t-H ' numerical simulations that the observed inhomogeneity of two-point correlations of galaxies 

can be accounted for by these models. Also, the natural appearance of bulk flows for an 
fT") . off-central observer is demonstrated. Finally, the inhomogeneity of the baryon content is 

discussed from the viewpoint of our inhomogeneous models. 

m : 

<N . 

y— I ■ §1. Introduction 

m ; 

' Under the assumption of spatial homogeneity, cosmological parameters have 

been determined using important observational results, such as the [magnitude m - 
redshift z] relation of type la Sm^JJifP 'Jj and the anisotropy of cosmic microwave 
Cu background radiation (CMB),eP'B)'LP'0)'eP and it has been found that models with 

a dominant cosmological constant best describe the SN data and flat models best 

', describe the CMB data. 

The present author I12P has considered models with a local void on scales of 
~ 200 Mpc in order to explain the existence of the large-scale bulk flows measured 
.,— i \ by Hudson et al. 0) and Willick,© and has shown that the zero A and 

small A models with flat space in the outer region are consistent with the data of 
type la supernovae (especially the recent data for z = 1.7E3 1 ) and with the CMB 
anisotropy. It thus appears that these models are competent for accounting for the 
cosmological observations in spite of the small realization probability of their large- 
scale void. Here the void means a low-density region with respect to the total matter, 
which consists mainly of the dark matter. 

The homogeneity of galactic distributions found in recent redshift surveys may 
support homogeneous models, but it does not necessarily rule out the inhomogeneous 
models with a local void, in which the galactic number densities in the inner and 
outer regions may be comparable owing to a larger suppression of galaxy formation 
in the outer region. 

In this paper we consider in §2 the process of galaxy formation and observa- 
tional aspects of inhomogeneous models from the viewpoint of hierarchical galaxy 
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formation theories. To demonstrate the feasibility of the inhomogeneous models, it is 
necessary to show how they can predict a nearly homogeneous distribution of galax- 
ies in spite of their void structure. Here we study the observed galactic distributions 
in connection with the existence of a local void and the two-point correlation that re- 
flects gravitationally the inhomogeneous structure of our models. In §3, we describe 
the models with a local void using numerical simulations, and in §4, we analyze the 
theoretical two-point correlations (£ and £ n ) in the inner and outer regions and 
show that the computed ratio £} is consistent with the observed ratio. In §5, we 
derive the bulk flows for an off-center observer. In §6 we discuss the inhomogeneity 
of the baryon content from the point of view provided by our models. Section 7 is 
dedicated to concluding remarks. 



§2. Galaxy formation and observational aspects 



Modern theories of galaxy formation led to a great deal of progress in the un- 
derstanding of galactic properties under the assumption of cold dark matter (CDM) 
in hierarchical clustering cosmologies. In these models, it is assumed that galaxies 
form when gas cools and condenses in dark matter halos that merge during their 
evolution. The key point in these theories is that the energy released from stars acts 
as a "negative feedback" on the gas and star formation. This idea was proposed by 
White and Rees® and developed by White and Frenk. c3) 

Recent works on galaxy formation, including the process of feedback, have em- 
ployed two important methods, a semi-analytic method and numerical simulations. 
The first is directly connected with the above two pioneering works and takes into 
account the feedback in detail, but has been applied only in simplified (spherically 
symmetric) situations. The merger history of a dark matter halo imsjpeen treated 
using either a statistical method (Kauffman et al.,HHP Cole et al.,^'™ Somerville 
and PrimackEj)) or with iV-body simulations (Kaufmann, Nusser and Steinmetz^^f 
Kauffmann, Colberg, Diaferio and White, c3 Benson et al.,c3) Somerville et al.c$). 
In the case of simulations, both the merger history (using iV-body simulations) and 
the feedback process and star formation (using gas dynamical simulations) have been. 
investigateiL(Cen and Ostriker,© Katz, Weinberg and Hernquist,© Frenk et al.,© 
Kay et al.E2P). This method has the advantage that no artificial symmetries need 
to be imposed, but it has the disadvantage that the treatment is not transparent 
and there are unrealistic limits on the dynamic range of resolved structures, though 
the resolution is being improved with the use of high-speed supercomputers. In the 
following, we first give an outline of recent treatments of galaxy formation on the 
basis of the semi-analytic method of Kauffmann et al., and next consider how we 
should treat galaxy formation in the present models with a local void. 

2.1. Galaxy formation in homogeneous models 

The semi-analytic models of Kauffmann et al.@) consist of a combination of 
./V-body simulations and an analytic approach to star formation and supernova feed- 
back. The ./V-body simulations are performed using the cluster normalization by 
determining <jg, the dispersion of density perturbations, within 8/i -1 Mpc spheres in 
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two CDM models: SCDM (rmodel) and ylCDM, with 

(J?o, A 0) h, o- 8 , r, f b ) = (1.0, 0.0, 0.5, 0.6, 0.21, 0.1), (0.3, 0.7, 0.7, 0.7, 0.9, 0.21, 0.15), 

respectively, where r is the shape parameter of the power spectrum and /& is the 
baryon factor (= Q^/Qq). 

The merging history of dark matter halos is constructed using the halo catalogues 
that consist of halos larger than the least stable system (containing 10 particles). In 
each halo, there is a central particle, representing the central galaxy onto which gas 
in a halo falls and where stars form. The central galaxy in a halo is the central galaxy 
in its most massive progenitor, and the central galaxies in less massive progenitors 
are satellites in the halo. These galaxies also grow from small galaxies to larger 
galaxies through merging. 

The physical properties of galaxies are determined by gas cooling, star formation, 
supernova feedback, dust extinction, etc. The most important processes among them 
are star formation and supernova feedback. Here, a star formation rate of the form 

dM^/olt (= M*) = aM co i d /t dyn (2-1) 

is assumed, where a is a free parameter, i<j yn is the dynamical time of the galaxy, 
and M* and M co id are the total masses of stars and cold gas in the halo, respectively. 

The energy ejected by a supernova explosion into the interstellar medium reheats 
the cold gas to the virial temperature. The reheating rate is expressed as 

dM reheat /dt = e^M,rj sn E sn /V c 2 , (2-2) 

where r/ sn is the number of supernovae per solar mass of stars (= 5 x 1O" 3 /M ), 
E sn is the kinetic energy of the ejection from each supernova (= 10 51 erg), e is a free 
parameter representing the fraction of this energy used to reheat cold gas, and V c is 
the circular velocity of galaxies. 

The evolution of distributions of dark matter halos and galaxies within them is 



determined using Eqs. (2-1) and (|2-2| ), when we specify the free parameters a and 
e, and the present distributions can be compared with observations in the form of 
luminosity functions, Tully-Fisher relations and two-point correlations. In order to 
choose the best values of a and e, a normalization condition is imposed: a fiducial 
reference galaxy (which is defined as a central galaxy with V c = 220 km s _1 ) is 
assumed to satisfy the /-band Tully-Fisher relation. Then the set (a, e) is given 
as (0.07,0.15) and (0.1,0.03) for the SCDM (rCDM) and ylCDM models, respec- 
tively. These values imply that the SCDM model has smaller star formation and 
larger feedback than the ACBM model. For these sets, both the SCDM and ACBM 
models have luminosity functions, Tully-Fisher relations and two-point correlations 
consistent with the observed ones. In addition to the feedback due to supernova 
explosions, photpipmzatipn and heating due to ultraviolet radiation also contribute 
to the feedback. H)©'© 

2.2. Galaxy formation in models with a local void 

In these models the total matter (consisting mainly of dark matter) has different 
uniform densities in the inner region (r < r&) and the outer region (r > r&), and, sim- 
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ilarly the luminous matter (consisting of galaxies) may also have different densities 
in the two regions. If the star formation and the feedback effect in the two regions 
act in a cooperative manner, however, the distributions of the luminous matter can 
be nearly uniform throughout the two regions. 

Let us assume for instance that the outer region is represented by the Einstein-de 
Sitter model and the inner region is represented by an open model with (Qq, Xq) = 
(0.3,0). Then we can have (0.07,0.15) as the set (a, e) in the outer region. Because 
the model in the inner region is similar to the low-density model with ({2q,Xo) = 
(0.3,0.7), we have (a, e) ~ (0.1,0.03) in the central part of the inner region, say 
r < r- m (< r b ). 

In the transient region satisfying n n < r < r&, the values of the set (a, e) change 
gradually from the value (0.07, 0.15) in the outer region to the value (0.1, 0.03) in the 
inner region. The strong photoionization by UV radiation from the outer region to 
the inner region may play an effective role in preventing small galaxies from forming 
and hence growing through merging in this region. For this reason, the number 
density of forming galaxies is small there, compared with the densities in the outer 
and inner regions. 

Galaxies formed in the transient region enter the outer region thereafter, colliding 
with other galaxies there, since the inner expanding velocities are larger than the 
outer velocities. At the present epoch, the boundary (r = r- m ) of the inner central 
region reaches the boundary r = of the outer region, so that the number densities 
of galaxies in the outer and inner regions may be nearly equal, as shown in Fig. 1. At 
present, however, we cannot determine the position r = rj n accurately, since studies 
on the dependence of rates of star formation, supernova explosion and UV radiation 
on the background models have not yet been established. Ej) Depending on these 
rates and the position of the central part (r < r; n ), the various large-scale structures 
may appear around the boundary (r = r&). It is found in simple simulations that, if 
radial distances for r- m < r < are about 1/4 of those satisfying < r < r& on the 
comoving scale, the two regions become smoothly connected later. 

at the early stage when galaxies are forming 
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i no feedback 
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Fig. 1. The schematic evolution of galactic density p g around the boundary rb- The dashed curves 
correspond to the case of no feedback. 

Here we have the important prediction that around the boundary there must be 
many active galaxies that are born as a result of the frequent collisions of galaxies, 
even if the distribution of galaxies is smooth there owing to strong feedback. 
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It is important to investigate the evolution of luminous matter using numerical 
simulations, taking into account the above described complicated processes including 
cooling, star formation, supernova feedback and heating by UV radiation, but such 
investigations are far beyond the scope of this paper. At present, we cannot show 
quantitatively how luminous matter becomes nearly homogeneous, but we are able 
to describe a possible process leading to the approximate homogeneity. 

2.3. Observed galactic distributions 

Here we consider the present observational situation regarding the galactic dis- 
tribution from the viewpoint of our inhomogeneous models. 

2.3.1. Dependence on the observer's position 

Several redshift surveys have suggested the existence of a low density xegkm or a 
local void with radius ~ 300/t Mpc (H = WOh km sec" 1 Mpc-^.SSSfl) The 
latest surveys (2dF and SDSS) show that the galactic distributions within a radius 
~ 600/i Mpc are homogeneous on the whole, though they include various large-scale 
structures. Thus the observed distributions seem to be consistent with homogeneous 
cosmological models. However, inhomogeneous models with a local void also may be 
compatible with these surveys for the following reasons: 

(1) The number densities of galaxies in the inner low-density region and the outer 
high-density region may be similar in these models also, as discussed in the previous 
subsection. 

(2) Some of observed large-scale structures, which are regarded as independent, may 
represent correlated structures near the boundary. For an observer in the center, 
these correlated structures have an equal redshift and can easily be distinguished. 
For us (off-center observers), however, they have direction-dependent redshifts and 
can be confused with other uncorrelated structures. When we superpose the redshift 
histograms in various directions to the north and south, off-central observers may 
see the structures there as independent at different distances. 

2.3.2. Active galaxies around the boundary 

From a recent study of spectral types of galaxies in the 2dF survey (Madgwick 
et al.t^iP), it was found that there are many active galaxies in the redshift region 
corresponding to the boundary r = r^. These galaxies may have been born as a 
result of the above-mentioned collisions of galaxies around the boundary. Similar 
analyses of the spectral types of galaxies are necessary for the SDSS survey to clarify 
the astrophysical situation around the boundary. 

2.3.3. Two-point correlation 

The apparent distribution of luminous matter may not be useful in the study of 
the background model consisting mainly of dark matter, because of the very compli- 
cated processes of formation and evolution of galaxies. The two-point correlations 
of galaxies, on the other hand, depend mainly on the background model, at least for 
less luminous galaxies in the homogeneous regions. Therefore it may be important 
to use these correlations to study the distribution of dark matter also. 

Recently, large-scale redshift surveys have provided valuable information about 
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the spatial dependence of the two-point correlations of galaxies in the region z < 
0.1. Using the NGP + SGP data from the 2dF surveys, Norberg et al.e2P derived 
the correlation length ro for samples with the various parameters ranges (cf. Ta- 
ble 1 and Fig. 3 in their paper). The observed correlation length increases with 
the absolute magnitude as r = 4.14 ± 0.64, 4.43 ± 0.45, 9.38 ± 1.48 for 

(M bJ - 51og 10 M) = (-18.0 18.5, 0.010 - 0.080), (-18.5 19.0,0.013 - 

0.104), • • • , (-21.5 - -22.5,0.059 - 0.28), respectively, where M bJ is the 6-band ab- 
solute magnitude. The SSRS2, EPS and Stromlo data also exhibit a trend similar 
to that found in the 2dF survey, as shown in Fig. 3 in Norberg et al.'s paper. c2P 
The luminosity dependence of vq in these samples seems to be consistent with that 
discussed in the hierarchical galaxy formation theories. 

According to a recent theory of galactic distribution (Benson et al.EP), ro is an 
increasing function of —M b j for — (M b j — 51og 10 h) > 21, while it is constant or a 
slightly decreasing function of —M b j for — (M b j — 5 log 10 h) < 21. The observed value 
of r is an increasing function of —M b j even for — (M b j — 51og 10 h) < 21. Quantita- 
tively, it increases by a factor of about 1.5 over the interval (M b j — 51og 10 h) = —18 

21. The two-point correlation £ cx rj (7 ~ 1.7) changes by a factor of 2.0 over 

this interval. We can interpret this change to represent the redshift dependence (or 
the spatial inhomogeneity) of ro or £ in the region z = 0.010 - 0.280. The boundary 
(z ~ 0.07) of the local void we consider is included in this region. 

From the data of the SDSS survey, Zehavi et al. derived the correlation length 
r for three samples, obtaining r = 7.42 ± 0.33, 6.28 ± 0.77 and 4.72 ± 0.44 for 
(M*, z) = (-23.0 - -21.5, 0.100 - 0.174), (-21.5 - -20.5, 0.052 - 0.097), and (-20.0 

18.5,0.027 - 0.051), respectively (cf. Table 2 and Fig. 16 in their paper.) Their 

result also reveals a change in ro by a factor of about 1.33 over the interval M* = 

—21.5 18.5. This change can be interpreted as the redshift dependence of ro in 

the interval z = 0.027 - 0.097. Our boundary (z ~ 0.07) is also included in this 
interval. 

In order to explain these spatial changes in ro and £, we study inhomogeneous 
models with a local void in the following sections using simple simulations. 

§3. Numerical inhomogeneous models 

In previous papers, I used models with a local void consisting of homogeneous 
inner and outer regions with a singular shell or an intermediate self-similar region. In 
order to study the evolution of nonlinear perturbations and their two-point correla- 
tions in these regions, I derived numerical inhomogeneous models using the method 
of iV-body simulations, by considering a spherical low-density region in the back- 
ground homogeneous models. These models are presented in this section. The back- 
ground model parameters are expressed as Hq, ao, f?o and Ao, representing the Hub- 
ble constant, the scale factor, the density parameter and the cosmological-constant 
parameter at the present epoch t = to, respectively, and the spatial curvature is 
K = OqHq{Qq + Ao — 1). The proper radial distance R at to is related to the ra- 
dial coordinate r by R = a^r. As the background models, we consider the following 
three cases: (J?o, Ao, h) = (1) Einstein-de Sitter model (1.0, 0.0, 0.5), (2) open model 
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(0.6,0.0,0.6) and (3) flat nonzero-yl model (0.6,0.4,0.6). 

The iV-body simulations were performed using the tree-code used by Suto and 
Suginohara, in which a periodic condition is imposed, and the initial perturbed state 
is determined using COSMICS. Here N = 2.1 x 10 6 , and the particle mass M is 
M = 5.7, 2.9 and 2.9 x 10 13 M Q for the above three cases (1), (2) and (3), respectively. 
The softening radius is 1 Mpc for all cases. The periodic condition is given at 



\x 2 \ = |x 3 | = r„, where R v = aor p = 300/h Mpc and r 



< p , VV11C1C ±Lp 

The initial conditions are set at t 



[(^1)2 + (a .2 )2 + ( a .3)2]l/2 i 



ti at which z = Zi (= 15.5), in the form 
of displacements 5x k and velocities v k of N particles (k = 1,2,3). After the initial 
conditions in the homogeneous case are given, a low-density region is introduced at 
the initial epoch by changing the particle positions x k in the inner region (r < r&) 
and the intermediate region (r^ < r < r^i) as 



x k — (x fc )h om x 



(1 + d) 



1 + d 



l 

rbi 



1 



for 
for 



r < r b , 

r b <r < rbi, 



(3-1) 



and x k = (x fc )hom for r > Hi (see Fig. 2). Here (2; fc )hom represents the particle 
positions in the homogeneous case and d is a constant expansion factor adjusted so 
as to give the expected density parameter and average expansion rate in the inner 
region. The intermediate region was introduced to make smooth the change in the 
density and velocity of particles near the boundary. Here, we treat mainly the case 
Rb = a^rb = 180/ h Mpc and Rbi = ao^bi = 210 /h Mpc, and consider also the case 
Rb = 120 /h Mpc and Rbi = 140/ h Mpc for comparison. 
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Fig. 2. Two regions and the boundary. The solid and dotted curves denote the surfaces with r = rt 
and r = r&i, respectively. 

From the viewpoint of our models with a local void, the above model parameters 
describeJhe_outer region (r > r&) and correspond to (^o 11 , Ao 11 , h 11 ) in the previous 
papers, Hy^'Ej^ where = lOO/i 11 km s _1 Mpc" 1 . From the simulations, the pa- 
rameters (f^o 1 , h l ) in the inner region are derived as their average values in the 
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expression 



#o^<+( 5> W /r (i) )/iV. 
iei* 



(3-2) 



along with Aq = \Aj (Hq) 2 , where the summation is taken over the set of particles 
included in the region r\ < r < r b (n ~ 0.2r b ), so as to avoid the disturbances 
near the origin and the boundary. Here, N* is the particle number in this region, 
and and are the radial velocities and radii of the z-th particle, where = 



k=i 



V?-\X 



(E 

For the particles that were initially in the inner region, the observed positions (in 
redshift space) are different from those represented by the background coordinates 
(in real space). To take the average velocity (in the inner region) into account, we 
define other coordinates (average comoving coordinates) x k (k = 1,2,3) as 



x k = x k x 



1 + 



Hi 



— J / 

r b r bl 1 1 



1 

m 



for 
for 



r < r b , 
r b <r < r bl , 



(3-3) 



and x k = x k for r > r b i ■ This coordinate system constitutes both an outer comoving 
system and an inner comoving system with respect to the mean motion. Using the 
volume V* = |7rao(rj) — rf)(H\/ Hq 1 ) 3 in these coordinates, the density parameter in 
the inner region is defined as 



nk = 



3ttG MN* 



(3-4) 



These values in the above three cases are listed in Table I together with d and R b - 
The values of d were chosen so as to obtain fi\ ~ 0.3 for three sets of (Oq 1 , Aq 1 ). By 
comparing the case R b h = 180 and the case R b h = 120, we find that the choice of 
the position of the boundary is not crucial. 

Table I. Inner model parameter values determined statistically for given outer model parameters 
(fla, Ao, h) — {fio 11 , Ao 11 , h u ). Here d is a constant expansion factor and Rb is the radial distance 
between the center C and the boundary, fib is the present baryon density parameter (oc J?o)- 





d 


R b h Mpc 


fio 1 Ao 1 


h l 




1.0 0.0 0.50 


0.022 


180 


0.38 0.00 


0.576 


1.98 


1.0 0.0 0.50 


0.022 


120 


0.39 0.00 


0.574 


1.98 


1.0 0.0 0.50 


0.030 


180 


0.32 0.00 


0.594 


2.21 


0.6 0.0 0.60 


0.040 


180 


0.29 0.00 


0.693 


1.55 


0.6 0.4 0.60 


0.040 


180 


0.29 0.29 


0.694 


1.55 



The distribution of dark matter particles at present epoch is displayed in Fig. 3 
as an example for the model parameter set (1.0,0.0,0.5), where the particles in the 
range —3 Mpc < aox 3 < 3 Mpc are plotted. The corresponding distributions of 
particles in the x k coordinates are displayed in Fig. 4. It is found that at the present 
epoch, the (dark matter) particles in the inner region near the boundary seem to 
have been mixed with those in the outer region, so that their distributions are quite 
complicated in the outer region near the boundary. 
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Fig. 3. The distribution of dark matter particles in the (aox 1 , aox 2 ) plane at the present epoch in 
the case {Ol 1 , Xl 1 ) = (1,0) and R b h = 180 Mpc. The length of the region displayed is 600/ h 
Mpc, and its width is 6/h Mpc. 
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Fig. 4. The distribution of dark matter particles in the (aox 1 , aox 2 ) plane at the present epoch in 
the case (fij 1 , Aj, 1 ) = (1,0) and R b h = 180 Mpc. The length of the region displayed is 6OO//1 
Mpc, and is width is 6/h Mpc. 
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§4. Theoretical two-point correlations 



In this section the two-point correlations (£) are derived from the numerical 
models given in the previous section. Their values in the inner region, with r < r^, 
and the outer region, with r p > r » r&i, have simple behavior, but they are 
complicated in the intermediate region, r& < r < r&i, because uncorrelated particles 
are mixed there. Here the correlations (£j and £ G ) in the inner and outer regions, 
respectively, were calculated using the particle positions in the average comoving 
coordinates (x k ) for various cases with R^h = 180 Mpc, which were treated in the 
previous sections. They are assumed to take the form £ = (ro/i?) 7 with 7 = 1.7 to 
derive the correlation lengths ro. The results are given in Table II. 

Table II. Correlation lengths and ratios of two-point correlations in the two regions of the models 
with R b h = 180 Mpc. Here, h = ft" 



(2o Ll 


AO 


h l /h u 




rift 




n l l ti l /{n l ti) 


1.0 


0.0 


1.152 


5.2 


4.2 


1.5 


2.3 


1.0 


0.0 


1.188 


5.6 


3.8 


2.1 


2.6 


0.6 


0.0 


1.155 


6.5 


5.0 


1.6 


1.8 


0.6 


0.4 


1.156 


5.4 


3.9 


1.8 


1.8 



It is found from this result that the ratio £ /£i is about 2.0 in the case l?o n = 
1.0, An 11 = 0.0, h l /h 11 = 1.188, and in other cases, it is somewhat smaller. Thus the 
situation in which the ratio of correlations for dark matter particles is larger than 
1.5 is common to these models with a local void. Here we compare this ratio with 
the observed ratio for galaxies. In general this is difficult, because the clustering of 
galaxies depends not only on the background model, but also on their luminosity 
through the bias effect. For galaxies with comparatively low luminosity, however, 
the bias effect on the clustering seems to be small, and the ratio for dark matter 

may be equal to that for galaxies with M^j = —18 21. Thus it can be concluded 

that the theoretical ratio (£ /£i) can account for the change in the observed ratio for 
galaxies over the range M^j = —18 21. 

Now let us consider the power spectra in the two regions. If the spatial scales 
are much less than r& or much larger than rt,, we can treat the power spectrum in 
each region as that of a homogeneous model with the fitting formula (i = I and II) 

"ln(l + 2.3% 



Pi(k) oc k 



2.3% 



[1 + 3.8% + (lG.lqiY + (5.46g0 + (6.71%) 4 ] 



-1/2 



(4-1) 



where qi = k(h l ) and the shape parameter is 1^ = Q^h 1 exp[— f] l B (l+y 2h l / Qq)], 

depending on the baryon density. EH In the outer region we can derive O^b} 1 from 
the measurements of CMB anisotropies, but at present we have only the rough value 
fi^h^ 1 = 0.3 - 0.5, which is obtained as those with weak prior from Boomerang, 
DASI and MAXIMA experiments. B)»0>»B In the inner region also we can roughly 
determine ft^h 1 from the treatment of Peacock and Dodds,!^ in which they derived 
the reconstructed linear data for the density contrasts of galaxies and clusters. Their 
data for the case k/h < 0.02 Mpc -1 correspond to our outer region and have com- 
paratively large uncertainty. Therefore, by assigning larger weights to the data for 
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k/h > 0.03 Mpc -1 , we estimate fi\h} as fi\h} = 0.2 - 0.25. Thus we obtain the 
observational ratio J7q I /i ii /(1?q/i i ) = 1.0 - 2.5 at present. 

In our models, by comparison, we have the ratios listed in the last column of 
Table II in the cases (1), (2) and (3). It is seen that our models correspond to the 
case with comparatively large ratios. 



§5. Bulk flows in the inner region 



If we observe from the center C, the average velocity field would be isotropic, 
and the directions of the velocities of particles would be radial at each point in the 
inner region. For an off-center observer O (being a realistic observer), however, the 
average velocity field is anisotropic, and the non-radial component is found to be 
significant, when each velocity(?/) in the inner region is divided into the component 
(v r ) in the radial direction O — > A and the component (v p ) in the direction C — > O 
as 

v = v r I + v p n. 

Here I and n are the unit vectors in the directions of O — > A and C 
(see Fig. 5), and we have 



(5T) 
O, respectively 



and 



vl — (vn)(ln) 



vn — (vl)(nl) 



1 



(In? 



{Inf 



Then the mean physical value V p for v p is defined by 



«o 



E 

iel* 



W) 



/A*, 



(5-2) 



(5-3) 



(5-4) 



where i is the particle index and A* is the particle number in the set I*. This value 
depends on the model parameters and is proportional to the distance R co (= a-or co ) 
from C to O. By choosing O in Ao different directions with equal R co , we can derive 
the Ao independent values of V p . Using them, we obtained an average value of V p 
and the dispersion (a p ) for Ao = 6. Their values are listed in Table III for models 
with Rj, = 180/Zi Mpc and R co = 50 or 60 /h Mpc, where h = h 11 . The values for 
Vbulk = (Hq — H^)R co = 100(/i : — /i n )i? co also are listed, and it is found that V p is 
equal to Vbuik- In the inner region, the average velocity field of galaxies is equal to 
that of dark matter, and it therefore appears that this velocity corresponds to the 
bulk flow of clusters with the velocity ~ 700km s _1 , which was measured by Hudson 
et slx3 and Willick.EHP If the measured velocity is smaller, T4mik can be adjusted to 
it by taking smaller R co , because ^uik is proportional to i? co . 



§6. Baryon content 



Let us now consider the baryon content in our models with a void. Here we 
assume that all inhomogeneities evolved gravitationally from perturbations with very 
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Table III. Bulk velocities and their dispersions 
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center observer boundary 

Fig. 5. Velocity vectors in the inner region. 

small amplitudes at the stage of nuclear synthesis and that dark matter and baryons 
are well mixed, so that the ratio of the baryon density to the matter density is equal 
everywhere. Then in the case (J?q, i?^ 1 ) = (0.3, 1.0) and h 11 = h 1 x 0.8 in the inner (I) 
and outer (II) regions, we have = 1.0/0.3, so that ^"(/i 11 ) 2 = 2.1 x Ol(h 1 ) 2 

and Of{h 11 ) 2 = 2.1 x Q\{h 1 ) 2 . The ratios Q l ^{h n f / Q\(h l f for the models we used 
in our simulations are given in the last column of Table I. It is found from this table 
that in models with = 1.0, the ratio can be about 2, and in models = 0.6, 
the ratios are somewhat smaller. 

From observational studies of light elements, it has been found that in the re- 
mote past, we have Q^h? = 0.025, which can be derived from the analysis of deu- 
terium abundance in QSOs for z ~ 2. e^>e3>>e9) Moreover, recent studies of the CMB 
anisotropy suggest similar baryon content, Q^h 2 = 0.02 - 0.03.B''EP Using these 
observational results in the outer region and the above ratios in the two regions, 
we obtain f?bh 2 = 0.010 - 0.015 in the inner region. As long as they are produced 
adiabatically, the baryon/photon ratio (77) is equal and constant (~ 7 x 10 -10 ) every- 
where, and therefore, at the stage of primordial nucleosynthesis, the photon number 
densities also were inhomogeneous. At the present stage, however, they have been 
equalized in both regions owing to the free propagation oLphotons. 

If the so-called crisis of big-bang nucleosynthesis e^'S is actual and 77 is really 
inhomogeneous (i.e. r\ ~ 3ji 10 - 10 in the inner region) in connection with the 
primordial abundance of Li 7 ,Ej)'llZP it may be necessary in order to avoid this crisis, 
to invoke that the void was produced as a special large-scale structure including 
non-adiabatic perturbations. 



r 
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§7. Concluding remarks 



In this paper, we first considered galaxy formation in inhomogeneous models 
with a local void. The observed distribution of galaxies seems to be homogeneous, 
but we discussed the possibility that this homogeneity may be due to the complicated 
feedback processes of galaxy formation and evolution in an inhomogeneous distri- 
bution of dark matter halos. In order to clarify observationally the existence of the 
boundary of the void, it is therefore necessary to investigate galactic distributions 
in greater detail, such as the distribution of active (colliding) galaxies and two-point 
correlations. 

Next, we derived numerical inhomogeneous models consisting of dark- matter 
particles with a spherical low-density (inner) region and studied the two-point cor- 
relations in the inner and outer regions. It was found that they can be made consis- 
tent with the redshift dependence of observed correlations of low-luminosity galaxies. 
Moreover we derived the bulk flow found by an off-center observer, which may corre- 
spond to observed large-scale bulk flow. Finally, we studied the baryon content and 



discussed the possibilitvo: 
In previous papers i3' 



avoiding the crisis of light elements. 
US' 113 we calculated the luminosity and angular-diameter 
distances to derive the [m, z] relation, and found that the distances in our models with 
a local void are similar to those in /1-dominant homogeneous models in the redshift 
interval < z < £ m (~ 1). Thus in the observations using mainly the distances, such 
as the number-counts of galaxies and clusters, the theoretical relations N{m) in the 
above two different types of models are similar in the interval < z < z m . Even if 
the observations of N(m) rule out the SCDM model (A = 0),©'H' our models with 
a local void and A = may be consistent with the observations, like the yl-dominant 
homogeneous models. 

It is well known that we are in the low-density region (J?q ~ 0.3), but the 
total matter density and baryon content in the remote region {z > 1) have not 
yet been established clearly through observation. From the observations of the CMB 
anisotropies, we can constrain the model parameters and baryon content. The recent 
results of Boomerang, DASI and MAXIMA experiments, i^' %%> however, seem to 
differ and to be partially inconsistent with respect to the density parameter and the 
baryon content. Therefore a precise measurement by the MAP satellite is needed 
to determine these values more accurately. If the additional data from SNIa with 
z > 1.5 are also obtained from the SNAP satellite and large telescopes, like Subaru 
etc., the model parameters in the remote region could be determined with much 
better precision. 
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